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There are important differences between the 
latitudinal structure of trends for the second half 
of the 20th century and for the 21st century 
(2000-2014) (Fig. 3). For example, the Arctic lati- 
tudes have shown strong warming trends both 
over the land and ocean since 2000, but during 
the latter half of the 20th century, the ocean 
trends in this area are near zero. The longer-term 
50-year trend has more consistency in the rates 
of warming across all latitudes, and this is even 
more evident over the full period of record back 
to 1880 (fig. S1). There is a distinct Northern 
Hemisphere mid-latitude cooling in LST during 
the 21st century, which is also showing up in 
cooling of the cold extremes, as reported for the 
extreme minimum temperatures in this zone in 
(27). Atmospheric teleconnections and regional 
forcings could be relevant in understanding these 
short time-scale zonal trends. It is evident that in 
most latitude bands, the global trends in the past 
15 years are comparable with trends in the pre- 
ceding 50 years. 

Last, we considered the impact of larger warm- 
ing rates in high latitudes (24) on the overall 
global trend. To estimate the magnitude of the 
additional warming, we applied large-area inter- 
polation over the poles using the limited observa- 
tional data available. Results indicate that, indeed, 
additional global warming of a few hundredths of 
a degree Celsius per decade over the 21st century 
is evident (Fig. 1), providing further evidence 
against the notion of a recent warming “hiatus” 
(supplementary materials). 

Newly corrected and updated global surface 
temperature data from NOAA’s NCEI do not 
support the notion of a global warming “hiatus.” 
As shown in Fig. 1, there is no discernable (sta- 
tistical or otherwise) decrease in the rate of 
warming between the second half of the 20th 
century and the first 15 years of the 21st century. 
Our new analysis now shows that the trend over 
the period 1950-1999, a time widely agreed as 
having significant anthropogenic global warm- 
ing (D, is 0.113°C decade”, which is virtually 
indistinguishable from the trend over the period 
2000-2014 (0.116°C decade’). Even starting a 
trend calculation with 1998, the extremely warm 
El Niño year that is often used as the beginning 
of the “hiatus,” our global temperature trend 
(1998-2014) is 0.106°C decade ‘—and we know 
that is an underestimate because of incomplete 
coverage over the Arctic. Indeed, according to our 
new analysis, the IPCCs (7) statement of 2 years 
ago—that the global surface temperature “has 
shown a much smaller increasing linear trend 
over the past 15 years than over the past 30 to 
60 years”—is no longer valid. 
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BRAIN CIRCUITS 


A parvalbumin-positive excitatory 
visual pathway to trigger fear 


responses in mice 


Congping Shang,”” Zhihui Liu,’ Zijun Chen,”” Yingchao Shi,” Qian Wang," 


Su Liu,’ Dapeng Li,’ Peng Cao™* 


The fear responses to environmental threats play a fundamental role in survival. Little is 
known about the neural circuits specifically processing threat-relevant sensory information 
in the mammalian brain. We identified parvalbumin-positive (PV*) excitatory projection 
neurons in mouse superior colliculus (SC) as a key neuronal subtype for detecting looming 
objects and triggering fear responses. These neurons, distributed predominantly in the 
superficial SC, divergently projected to different brain areas, including the parabigeminal 
nucleus (PBGN), an intermediate station leading to the amygdala. Activation of the PV* 
SC-PBGN pathway triggered fear responses, induced conditioned aversion, and caused 
depression-related behaviors. Approximately 20% of mice subjected to the fear- 
conditioning paradigm developed a generalized fear memory. 


nvironmental threats are detected by dif- 
ferent sensory organs projecting to central 
brain areas to trigger fear responses (J, 2). 
The superior colliculus (SC) is a retinal re- 
cipient structure (3, 4) composed of different 
neuronal subtypes (5, 6), including parvalbumin- 
positive (PV*), somatostatin-positive (SST), and 
vasoactive intestinal peptide-positive (VIP*) neu- 
rons (Fig. 1A and fig. S1). In addition to mediat- 
ing orienting responses (7), the SC contributes 
to avoidance and defense-like behaviors (8-11). 
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With an optogenetic approach (12-14), we found 
that activation of neurons expressing channel- 
rhodopsin-2 (ChR2) in mouse SC triggered freez- 
ing that lasted 52.8 + 5.3 s (n = 5 mice) (movie 
S1). This prompted us to systematically iden- 
tify the key neuronal subtypes underlying this 
behavior. 

By crossing Ai32 (15) with different Cre lines 
(Fig. 1B) (16, 17), we expressed ChR2-enhanced 
yellow fluorescent protein (EYFP) in specific 
neuronal subtypes in the SC (Fig. 1C and fig. S1) 
and optogenetically elicited spikes in acute slices 
(Fig. 1D and fig. S1). Activation of SC PV* neurons, 
but not SST* or VIP* neurons, triggered impulsive 
escaping (1.18 + 0.09 s) followed by long-lasting 
freezing (46.4 + 2.8 s) (Fig. 1, E to G; fig. S1; and 
movie S2). To avoid activation of PV* retinal 
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ganglion cells (RGCs) (18) by ferrule light, we in- 
jected adeno-associated virus (AAV) expressing 
double-floxed ChR2-mCherry (12) into the SC of 
PV-ires-Cre mice, resulting in specific expres- 
sion of ChR2-mCherry in SC PV* neurons but 
not in PV* RGCs (Fig. 2A and fig. S2). The light 
triggered spikes from ChR2-mCherry-positive 
neurons in SC slices (Fig. 2B and fig. S2), elic- 
ited a similar stereotyped locomotor pattern 
(fig. S2 and movie S3), and increased the heart 
rate and plasma corticosterone levels that were 
not observed in mice with SC PV* neurons ex- 
pressing mCherry (Fig. 2, C and D). 

When facing threats, animals can either fight 
or flee. To test whether SC PV* neurons were 
involved in this behavioral dichotomy, we mea- 
sured the durations of light-induced escaping 
(E) and freezing (F) and calculated their ratio 


(E/F ratio). We conducted a series of tests span- 
ning 5 days (table S1). First, light stimulations 
with higher intensity or longer duration enhanced 
E/F ratios in the same male mice by prolonging 
escaping more strongly than freezing (Fig. 2, E 
and F). Second, light stimulations with higher 
frequency but similar total illumination time 
prolonged escaping and freezing proportionally 
(Fig. 2G). Third, both responses showed strong 
adaptation to repetitive light stimulations (every 
5 min), with no significant change in E/F ra- 
tios across each stimulation (Fig. 2H). Finally, 
the same light stimulations elicited longer es- 
caping and shorter freezing in female versus 
male mice, resulting in higher E/F ratios in fe- 
males (Fig. 21). The origin of these sexually 
dimorphic behaviors was further examined (sup- 
plementary text). 


We next characterized the morphological and 
physiological properties of SC PV* neurons. 
They were predominantly but not exclusively 
distributed in the superficial gray (SuG) layer of 
the SC (Fig. 3A and fig. S3). Whole-cell record- 
ing of tdTomato-expressing PV* neurons in SC 
slices from PV-ires-Cre; Ai9 mice (19) demon- 
strated that, distinct from V1 PV* interneur- 
ons with slow frequency adaptation (20), the 
SuG PV* neurons responded to depolarizing cur- 
rents in a faster adaptation mode (Figs. 3, C 
and D, and fig. S3). SuG PV* neurons labeled 
with neurobiotin had parallel dendrites extend- 
ing to the SC surface, presumably receiving 
inputs from RGCs (Fig. 3B and fig. S3). The 
postsynaptic currents from PV-negative neurons 
induced by optogenetic activation of PV* neu- 
rons expressing ChR2-mCherry were blocked by 


/ 2 ms 10Hz 
x < 
S S 
h ® 
8 g 
> > 
a 1mm 250 pm a 
B 2 ms 10Hz 
NN / Sy = m = = 
roter < Cay E 
Ai32 STOP § S 
a nH 
g £ 
BS = 
Cre lines SST, PV or VIP gene [ies ] C ce} 5 5 
o mm 250um 2 0.1 s|20 mV 
E 12- F 
PV-ires-Cre Light stimulation Before Stimulation (5 sec) During Stimulation (5 sec) After Stimulation (30 sec) During Stimulation (5 sec) Stimulation (5 sec) 
T SST-ires-Cre PV-ires-Cre Start/End ® 1.2-7 tek 1.475 
2 0.9- a Eng Ẹ 1.0-4 5124 
£ nd Start Start g E io 
= £0.84 Soo] 
S 3 & 08+ 
= 2 06+ 2 
£ 0.64 SST-ires-Cre a a z 064 
8 a 3 0.4- & 0.4 
5 End Start|| End È aa] ns. F aie 
3 Start End ion i 
2 10 cm 0 0 x 
7) — E i & wv 
n G After Stimulation (30 sec) 
| I | I | á - z 
-5 0 5 10 20 25 30 35 ae “a 2 ag 
Time relative to light stimulation initiation (sec) 30104 fy z 3.04 pu ------ 
E 0.08 5 eg 
2.04 
0.06 | $ 
a 3457 
x 0.044 g 
D & 1.04 
2 0.024 ihe Fost 
9} 9 i ~ Jioll o ii 
x AX ox o AX oX 
Se 8 SS 


Fig. 1. Neuronal subtypes in the SC to trigger fear responses. (A) Diagram 
of different neuronal subtypes in the SC. (B) Ai32 mice were crossed with 
different Cre lines. (C) Coronal micrographs showing ChR2-EYFP expressed in 
specific neuronal subtypes. DAPI, 4’,6-diamidino-2-phenylindole. (D) Light- 
induced spikes from ChR2-EYFP* neurons in acute slices. (E) Instanta- 
neous locomotion speed before, during, and after light stimulation. (Inset) 
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Locomotion trails from example mice. (F and G) Analyses of peak speed 
and travel distance during and after light stimulation. Data in (F) and (G) 
are means + SEM (error bars); numbers of mice are in bars. Statistical 
analysis is t test (***P < 0.001; n.s. P > 0.1). Dashed lines in (F) and (G) 
indicate the control levels measured from behaviors after SC SST* neuron 
activation. 
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p-(-}-2-amino-5-phosphonopentanoic acid (APV) 
and 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), 
but not by picrotoxin (Fig. 3, E and F, and fig. 
S3), suggesting that these neurons were gluta- 
matergic but did not release y-aminobutyric 
acid. The SC PV* neurons in intermediate and 
deep layers of the SC were examined (supple- 
mentary text). 

To test whether SC PV* neurons were involved 
in detecting collision threats in the visual field 
(21-23), we displayed a virtual soccer ball moving 
in controlled velocities and directions to the 
anaesthetized mice (Fig. 3J and fig. S4). The single- 
unit activity recorded with optrodes was quan- 
titatively identified (24, 25) as putative SC PV* 
neurons expressing ChR2-mCherry (Fig. 3, G to 
I). These putative PV* neurons (n = 9 cells) were 
strongly activated by the ball moving toward 
the animal but not by the motion in the other 


five directions (Fig. 3K). The response onset 
time before collision depended on the size and 
velocity of the ball (Fig. 3L) and was linearly 
correlated with the square root of the diameter/ 
velocity (Fig. 3M). The response peak was close 
to the time to collision and was independent 
of the size and velocity of the soccer ball (Fig. 
3N). In freely behaving mice, the escaping trig- 
gered by SC PV* neuron activation pointed to 
the side of the SC receiving light stimulation 
(movie S4). 

We then determined the circuit mechanism 
underlying the fear responses mediated by SC 
PV* neurons. By injecting AAV expressing double- 
floxed monomeric green fluorescent protein 
(mGFP) into the SC of PV-ires-Cre mice, we spe- 
cifically labeled SC PV* neurons (Fig. 4B and 
figs. S5 and S6) and observed axon terminals in 
the parabigeminal nucleus (PBGN), the pontine 


nucleus (Pn), and the dorsal lateral geniculate 
nucleus (DLGN) (Fig. 4, A and C). These pro- 
jections were confirmed by retrograde tracing 
with cholera toxin B with Alexa Fluor-594 (CTB- 
594). CTB-594 injection into the PBGN (Fig. 4D) 
retrogradely labeled SC neurons predominantly 
in the ipsilateral SuG layer (Fig. 4E and fig. S7). 
A considerable proportion of CTB-labeled SC 
neurons (SC-PBGN: 52 + 4.3%; SC-Pn: 31 + 4.5%; 
SC-DLGN: 33 + 3.8%, n = 3 mice) were positive 
for PV (Fig. 4F and fig. S7). 

To determine which of these parallel projec- 
tions (26) participated in the fear responses, we 
injected AAV expressing double-floxed ChR2- 
mCherry into the SC and implanted optic fibers 
in the PBGN or Pn (Fig. 4G) to locally stimulate 
the ChR2-mCherry-positive axon terminals. Ac- 
tivation of the PV* SC-PBGN pathway, but not the 
PV" SC-Pn pathway, triggered the stereotyped 
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Fig. 2. Specific activation of SC PV* neurons induced fear responses. 
(A) Specific expression of ChR2-mCherry in SC PV* neurons of PV-ires-Cre 
mice. (B) The light-pulse train triggered spikes (red) from ChR2-mCherry-— 
positive neurons and postsynaptic currents (black) from adjacent ChR2- 
mCherry—negative neurons. (C) Electrocardiographic traces and heart rate 
analyses from the anaesthetized mice before and after light stimula- 
tions. Ctrl, control. (D) Analyses of plasma corticosterone concentration in 
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M, male; F, female. 


response to light stimulation. (E to 1) Durations of escaping, freezing, and 
E/F ratios, were plotted as functions of stimulation intensity, duration, 
frequency, repetition, and sex in mice with SC PV* neurons expressing 
ChR2-mCherry. Data in (C) to (1) are means + SEM (error bars); numbers 
of mice are in bars. Statistical analysis is t test (***P < 0.001; **P < 0.01; 
n.s. P > 0.1). Dashed lines indicate the levels measured from control mice. 
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escaping-freezing locomotor pattern (Fig. 4, H to 
J; fig. S9, and movie S5). We examined whether 
PBGN projected to the amygdala by anterograde 
and retrograde tracings. Local injection of AAV- 
SynaptoTag (27) in the PBGN and its adjacent 


region strongly labeled axon terminals positive 
for synaptobrevin-2-EGFP in the central amyg- 
daloid nucleus (28, 29), whereas CTB-594 injection 
in the amygdala retrogradely labeled neurons in 
the PBGN (fig. S8). Finally, the relation between 


PV* SC-PBGN pathway activation and the affec- 
tive state of mice was explored (supplementary 
text and figs. S10 to S12). Taken together, these 
data revealed a PV" excitatory visual pathway to 
trigger stereotyped fear responses in mice. 
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Fig. 3. Morphological and physiological properties of SC PV* neurons. 
(A) Layer-specific distribution of SC PV* neurons. (B) Neurobiotin-labeled 
PV* neurons in the SC SuG layer and V1. (C) Spike firings of PV* neurons in 
the SC and V1 to depolarizing currents. (D) Analyses of spike number and 
spiking phase as a function of current intensity. n, number of cells. (E and F) 
Effects of CNQX (20 uM)/APV (50 uM) (E) and picrotoxin (50 uM) (F) on the 
postsynaptic currents (PSC) induced by light stimulation. n, number of cells. 
(G) Single-unit activity recorded from a putative SC PV* neuron triggered by 
light pulses (arrows, 1 ms at 10 Hz). (H) Raster plot showing the latency of 
light-evoked spikes relative to the light pulses (O ms). (1) Distributional plot 
(left) and example spikes (right) evoked by visual stimuli and light showing 
quantitative identification of PV-positive and PV-negative units based on the 
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Square root of diameter/velocity 


waveform correlation and energy of light-evoked spikes. (J) A virtual soccer 
ball flying toward the eye of an anaesthetized mouse. (K) Example single-unit 
traces from a putative SC PV* neuron in response to the soccer ball (20 cm 
in diameter) moving in six directions at 2 m/s. (L) Peristimulus time his- 
tograms of a PV* neuron to looming stimuli with controlled velocities (V) (2 
and 5 m/sec) and diameters (D) (20 and 40 cm). Arrows indicate response 
onset time. (M and N) Correlation analyses of response onset time (M) and 
response peak time (N) of SC PV* neurons and the square root of 
diameter/velocity of the looming ball. Data are means + SEM (error bars); 
numbers of cells or units are in graphs. Statistical analyses are t test and 
one-way analysis of variance (***P < 0.001; n.s. P > 0.1). R, correlation 
coefficient. 
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Our data lead to the following conclusions. 
First, the SC PV* neurons form a subcortical visual 
pathway that transmits threat-relevant visual 


A B 
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information to the amygdala to trigger fear re- 
sponses. These data, in alliance with earlier studies 
(18, 30, 31), suggest a “retina-SC-PBGN-amyegdala- 


l After Stimulation (30s) 

0.125 3.6- MSk 
CLAE Ry ns. £307 
3 S 
£ 0.08- 8 2.4- 

5 S 

2 0.06- v 1.8- 

F 2 

x 0.047 2 1.2- 

E kkk kkk g 

0.024 9 9 F 0.6- o r 

11 m m 10 0 ill oo E 10) 
SC_SC_PBG Pn SC SC PBG Pn 


Ctrl ChR2 Ctrl ChR2 


Fig. 4. PV* SC-PBGN pathway mediated fear responses. (A to C) Specific 
expression of mGFP in SC PV* neurons (B) of PV-ires-Cre mice resulted in 
labeling of their axon terminals in the PBGN, Pn, and DLGN [(A) and (C)]. 
(D to F) CTB-594 injected in the PBGN (D) retrogradely labeled cells in 
the SC (E), a large proportion of which were PV* (denoted by arrowheads) 
(F). DpG, deep gray layer; InW, intermediate white layer; Ing, intermediate 
gray layer; Op, optic nerve layer; Zo, zonal layer. (G) Diagrams showing the 
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hypothalamus” pathway for vision-induced fear 
responses. Second, the SC PV* neurons in the 
SuG layer are predominantly glutamatergic 
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projection neurons with spiking patterns dis- 
tinct from those of their counterparts in cor- 
tical regions. Thus, this finding broadens the 
concept of PV* neurons (32) and adds another 
perspective to understanding their functions. 
Third, the SC PV* neurons may belong to type-p 
looming detector, supporting the notion that 
mathematically defined computational units 
correspond to specific neuronal subtypes (33). 
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A Cas9-guide RNA complex 
preorganized for target 


DNA recognition 


Fuguo Jiang,’ Kaihong Zhou,” Linlin Ma,” Saskia Gressel,® Jennifer A. Doudna’”**?©* 


Bacterial adaptive immunity uses CRISPR (clustered regularly interspaced short 
palindromic repeats)—associated (Cas) proteins together with CRISPR transcripts for 
foreign DNA degradation. In type II CRISPR-Cas systems, activation of Cas9 endonuclease 
for DNA recognition upon guide RNA binding occurs by an unknown mechanism. Crystal 
structures of Cas9 bound to single-guide RNA reveal a conformation distinct from both 
the apo and DNA-bound states, in which the 10-nucleotide RNA “seed” sequence required 
for initial DNA interrogation is preordered in an A-form conformation. This segment of 
the guide RNA is essential for Cas9 to form a DNA recognition—competent structure 

that is poised to engage double-stranded DNA target sequences. We construe this as 
convergent evolution of a “seed” mechanism reminiscent of that used by Argonaute 


proteins during RNA interference in eukaryotes. 


RISPR-Cas proteins function in complex 

with mature CRISPR RNAs (crRNAs) to 
identify and cleave complementary target 
sequences in foreign nucleic acids (J). In 

type II CRISPR systems, the Cas9 enzyme 
cleaves DNA at sites defined by the 20-nucleotide 
(nt) guide segment within crRNAs, together with 
a trans-activating crRNA (tracrRNA) (2) that forms 
a crRNA:tracrRNA hybrid structure capable of 
Cas9 association (3). Once assembled on target 
DNA, the Cas9 HNH and RuvC nuclease domains 
cleave the double-stranded DNA (dsDNA) sequence 
within the strands that are complementary and 
noncomplementary to the guide RNA segment, 
respectively (3, 4) (Fig. 1A). By engineering a syn- 
thetic single-guide RNA (sgRNA) that fuses the 
crRNA and tracrRNA into a single transcript of 
80 to 100 nt (Fig. 1B), Cas9:sgRNA has been har- 
nessed as a two-component programmable system 
for genome engineering in various organisms (5, 6). 
The utility of Cas9 for both bacterial immunity 
and genome engineering applications relies on 
accurate DNA target selection. Target choice relies 
on base pairing between the DNA and the 20-nt 
guide RNA sequence, as well as the presence of a 
2- to 4-base pair (bp) protospacer adjacent motif 
(PAM) proximal to the target site (3, 4). The tar- 
get complementarity of a “seed” sequence within 
the guide segment of crRNAs is critical for DNA 
recognition and cleavage (7, 8). In type II CRISPR 
systems, Cas9 binds to targets by recognizing a 
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PAM and searching the adjacent DNA for com- 
plementarity to the 10- to 12-nt “seed” sequence 
at the 3’ end of the guide RNA segment (Fig. 1B) 
(3, 9-11). Crystal structures of Cas9 bound to 
sgRNA and a target DNA strand, with or without 
a partial PAM-containing nontarget strand, show 
the entire 20-nt guide RNA segment engaged in 
an A-form helical interaction with the target 
DNA strand (12, 13). How the “seed” region with- 
in the guide RNA specifies DNA binding has re- 
mained unknown. 

To determine how Cas9 assembles with and 
positions the guide RNA prior to substrate recog- 
nition, we solved the crystal structure of catalyt- 
ically active Streptococcus pyogenes Cas9 (SpyCas9) 
in complex with an 85-nt sgRNA at 2.9 A reso- 
lution (Fig. 1 and table S1). The overall structure 
of the Cas9-sgRNA binary complex, representing 
the pre-target-bound state of the enzyme, resem- 
bles the bilobed architecture of the target DNA- 
bound state, as observed in electron microscopic 
studies (14), with the guide segment of the sgRNA 
positioned in the central channel between the 
nuclease and helical recognition lobes (Fig. 1, C 
to E). This structural architecture and guide RNA 
organization is maintained in the crystal structure 
of a widely used nuclease-inactive version of Cas9 
(D10A/H840A, referred to as dCas9) in complex 
with sgRNA (fig. S1). 

Comparison of SpyCas9 crystal structures rep- 
resenting the protein alone and the RNA-bound 
and RNA-DNA-bound states of the enzyme reveals 
the nature of Cas9’s conformational flexibility dur- 
ing sgRNA binding and target DNA recognition 
(Fig. 2A and figs. S2 and S3). The helical recog- 
nition lobe undergoes substantial rearrangements 
upon sgRNA binding but before DNA association, 
especially in helical domain 3, which moves as 
arigid body by ~65 A into close proximity with 
the HNH domain (fig. S2D). Superposition of 
the Cas9-sgRNA pre-target-bound complex onto 
the target DNA-bound structures reveals further 
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